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Abstract 

Recent experimental and theoretical studies on the magnetization dynamics driven by an elec- 
tric current have uncovered a number of unprecedented rich dynamic phenomena. We predict an 
intrinsic chaotic dynamics that has not been previously anticipated. We explicitly show that the 
transition to chaotic dynamics occurs through a series of period doubling bifurcations. In chaotic 
regime, two dramatically different power spectra, one with a well-defined peak and the other with 
a broadly distributed noise, are identified and explained. 

PACS numbers: 75.75.+a, 85.75 .-d, 05.45.-a, 05.45.Pq 
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Among many well-studied non-linear oscillators driven by external forces, only a hand- 
ful oscillators have technological applications l[ 2]. The recently discovered current-driven 
magnetization oscillators with tunable microwave frequencies in spin valves are very desir- 
able for magnetic storage devices and for telecommunications. Up until now, theoretical and 
experimental studies of the oscillator have been carried out only for the simplest cases where 



the dynamics of the oscil 



ator is either in a self-sustained steady-state precessional motion 



(limit cycle) (j], 4], l], [(3, B, |j] or in synchronization with other oscillator(s) jsl, [l(3, 11, Q, l^ - 
Since the equation that governs the oscillator is highly non-linear, it would be fundamentally 
interesting to map out the full dynamics for experimentally relevant parameters. In partic- 
ular, a thorough study of chaotic dynamics will elucidate how the current-driven oscillator 
responds to an external perturbation. 

Previous study of the current-driven magnetization chaos or noise was based on the mi- 
cromagnetic simulation where the magnetization is not uniform due to strong magnetostatic 
interaction at the edge of the sample These chaotic dynamics highly depend on 

the shape and size of the sample, and thus it is not an intrinsic property of the current- 
driven oscillator. Here we consider a single-domain current-driven spin valve oscillator so 
that the undesired complication of the spatial variation of the magnetization is eliminated. 
The intrinsic dynamic property of the spin valve oscillator is investigated in the presense of 
an external periodic perturbation, for example, an AC current. By utilizing the Poincare 
map we have found the route to chaos to be via period doubling bifurcations. Positive 
Lyapunov exponents fl^l] and Sharkovskii ordering _17 1 are observed as the evidence of chaos. 
Furthermore, we show two dramatically different power spectra in chaotic regions, one with 
a well-defined peak and the other with broadly distributed noise. 

The modeled spin valve consists of a pinned layer whose magnetization is fixed along the 
positive x-axis and a single-domain free layer whose magnetization vector m is the subject 
of our calculation. The free layer experiences an effective field H e ff made of an external 
field, an anisotropy and a demagnetization field perpendicular to the layer (z-axis). We 
choose the direction of the magnetic field along the in-plane easy axis The dynamics of 
the magnetization on the free layer is determined by the modified Landau-Lifshitz-Gilbert 

n 

(LLG) equation. |3j 

dm dm 

— = -7m x H eff + am x — + 0,7m x (m x e x ) (1) 
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where aj is the amplitude of the spin torque which has the unit of the magnetic field (Oe). 
Since the magnitude of the magnetization |m| = 1 is a constant, there are only two indepen- 
dent variables, thus, chaotic dynamics are excluded for a constant aj. In fact, the solutions 
of Eq. (pD) have already been obtained The time-dependent solution is a limit cycle 

which can be analytically determined via the construction of a Melnikov integral (MI) {§]. 
It has been found that the stable limit cycle with a well defined frequency (which we call 
the natural frequency Uq) exists only for the current density larger than a critical current 
density. Two distinct limit cycles have been identified: an out-of-plane orbit and a nearly 
in-plane orbit around m = e x . 

To explicitly reveal the dynamics of the above spin valve oscillator under an external 
perturbation, we now consider a time-dependent current which adds an additional term in 
Eq. (pQ), a ac cos(ut)jmx (raxej. The simplest dynamic phase would be the synchronization, 
i.e., the spin valve oscillator is forced to oscillate in phase with the external frequency uj as 
long as ujq is sufficiently close to w. Indeed, the recent experiment Q and theory [3, Q] 
have clearly demonstrated the synchronization. The much richer dynamics, however, would 
be chaotic dynamics {^(J shown below. 

An analytical prediction of chaotic dynamics can be made using a MI along the separatrix 
orbit. The simple zeros of the MI indicate the occurrence of chaotic dynamics 1 201. In 

n 

Ref. [20] . the MI was carried out for a fixed magnetic field. Here we map out the phase 
diagram in a parameter space of both DC current and magnetic field. Implementing the 
MI approach, we obtain a boundary loop which encloses all the possible simple zeros in the 
parameter space. In Fig. [T], the phase diagram is plotted. Two solid lines are a portion of 
the complete boundary loop; if we extend the parameter space to a much larger range, the 
two lines will meet to form a closed area. Since there exists an arbitrary phase in calculating 
the MI, the exact solutions are not available. Therefore, a more rigorous numerical test is 
required to identify the true chaos bounded by the boundary lines. We show below that 
only the dark region in Fig. [U is chaos. 

Numerically, chaos is often indicated by the positive Lyapunov exponents. The Lyapunov 
exponent measures the exponential increase (or decrease) of an initial distance from two close 
trajectories, 

x ,• 1 i \\ Sm t\\ 

Xi = hm - In- 



where A, is the ith Lyapunov exponent and ||5m^|| is the distance between the trajectories of 
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Figure 1: Chaotic/non-chaotic phase diagram. Two lines are the boundary of chaos determined by 
Melnikov Integral. Chaos is possible only within the region bounded by these two lines. The dark 
region represents a positive value of the largest Lyapunov exponent, i.e., chaos. The parameters 
are chosen similar to a permalloy film: damping constant a = 0.02, anisotropy constant Hk = 0, 
demagnetization field 4ttM s = 8400(Oe) and 7 = 1.7 x 10~ 7 (Oe _1 s _1 ). The amplitude and the 
frequency of the ac current are a ac = 20(Oe) and u = 15(GHz). 



the ith orthogonal axis at time t. When A« > 0, the distance between two arbitrarily close 
points will increase exponentially along the z-axis at a large t For a bounded dynamic 
system, any positive Lyapunov exponent indicates chaos. We numerically computed the 
Lyapunov exponents and Fig. [2](a) shows an example of the Lyapunov spectra, where the 
largest exponent is shown in red solid line and the other two are in dashed and dotted lines. 
When cij is smaller than the critical value 230.03(Oe), all three Lyapunov exponents are 
non-positive and thus the dynamics are non-chaotic, see Fig. [2](b). When a,j becomes larger 
than 230.03(Oe), the largest Lyapunov exponents become positive, and thus chaos appears 
[Fig. [2](c)]. By sweeping the parameters of the DC current and the external field, we have 
mapped out the parameters that give arise at least one positive Lyapunov exponent, shown 
as the dark area in Fig. [H 

To understand how chaos develops, we applied the Poincare map to this system. In 
spin valve oscillators, it is difficult to choose a simple Poincare section because a simple 
section can hardly include both out-of-plane and in-plane orbits simultaneously. To avoid 
this difficulty, we record the trajectory points every time a local minimum of m x is reached. 
If we define the time intervals between two successive minima as a period, this is essentially 
a Poincare period map except that the period is not a constant. 

By using the minimum m x map, we are able to see the development of chaotic dynamics 
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Figure 2: (a) Lyapunov spectra for 226(Oe) < a, < 233(Oe); (b) the regular non-chaotic trajectory 
for a,j = 227.5(Oe) ; (c) chaotic trajectory for Oj = 232.5(Oe). The magnetic field is 200 Oe. All 
parameters are same as in Fig. 1. 

when one varies the parameters. For limit cycles where the parameters are within the non- 
chaotic white regions shown in Fig. 1, the magnetization follows a unique trajectory and the 
map point is just one unique point. When the parameters, e.g., the currents, are close to 
the boundaries between the white and dark regions of Fig. 1, two map points are seen; this 
is identified as the period doubling bifurcation. In this case, the magnetization orbit will 
return to the original orbit after two periods. If we choose the current even closer to the 
boundaries, a period-four orbit and, in general, a period-2 n orbit, appear. When n — ► oo, the 
magnetization dynamics become chaotic, i.e., the appearance of the dark regions in Fig. 1. 
In Fig. [3j we have shown an example of the bifurcation diagram where the transition from 
the synchronization to the bifurcation and to chaos is clearly demonstrated. Period-five and 
period-six chaos windows are observed between 252 to 253 Oe in the cascade. According to 
Sharkovskii ordering theorem, these observations imply the existence of chaos [17]. 



The period doubling bifurcation cascade shown in Fig. [3] is in fact similar to those in other 



nonlinear systems 



16j. To see if our bifurcation structure belongs to the same universal 



quadratic class, we evaluate the Feigenbaum ratio 



22j which is defined as 



a i+l — Q-i+2 

where a-i is the value of the parameter at the iih bifurcation point. In Table fl], We list an 
example of measurement of Feigenbaum ratios for our system. There is a second Feigenbaum 




Figure 3: Period doubling bifurcation cascade. Dashed lines indicate the critical values of the DC 
current where the bifurcation occurs. All parameters are the same as in Fig. 2. 



Table I: Feigenbaum ratio measurement. The second column is derived from Fig. 3; at is a critical 
value at the ith bifurcation point. 
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CLi 




Si 


1 


257.2500 


-1.9300 


4.289 


2 


255.3200 


-0.4500 


4.356 


3 


254.8700 


-0.1033 


4.472 


4 


254.7667 


-0.0231 


4.529 


5 


254.7436 


-0.0051 


4.636 


6 


254.7385 


-0.0011 
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254.7374 







ratio that denotes the rescaling ratio of the bifurcation forks. For example, from first two 
bifurcation in Fig.[3l we measured \a\ = \AM x i/ AM x2 \ = 2.68. Analytically, the Feigen- 
baum ratios are approaching to the universal constants 4.6692. . . for 6 and 2.5029. . . for 



| a | [22j, for a quadratic non-linear map[23j. We find that the measured ratios agree with 
the universal Feigenbaum numbers exceedingly well. This indicates that our system indeed 
belongs to a quadratic nonlinear class. 

By summing all Lyapunov components, we have verified that the volume contraction (A„ = 
Ai + A2 + A3) is always negative even though the Lyapunov exponent can be positive; the 
negative volume contraction indicates the dissipative nature of the spin valve oscillator 
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A highly interesting feature is the two distinct magnetization trajectories during the tran- 
sition to chaos. In the first case, the bifurcation occurs only at the out-of-plane orbits shown 
in Fig. 0](a). Although the bifurcation on this single orbit also leads to chaos because the 
largest calculated Lyapunov exponent is positive, the power spectrum displays a well defined 
peak. The peak position corresponds to the inverse of the average time for the magnetization 
to complete one loop (since the loop never closes, we define a one-loop when the trajectory 
returns to the point nearest to the starting point of the loop). Thus, the presence of the 
narrow peak indicates a quasi-periodic motion of the magnetization in chaotic dynamics. It 
would be erroneous if one automatically assumes the dynamics is synchronization when the 
experimental power spectrum is highly peaked. Synchronization refers to the phase-locking 
between the external and natural frequencies but the positive Lyapunov exponent excludes 
the possibility of the phase-locking. The second chaotic motion involves both out-of-plane 
and in-plane orbits, as seen in Fig. IU(c). In this case, the power spectra display typical 
noise, i.e., broadly distributed spectra. The magnetization jumps between out-of-plane and 
in-plane orbits are completely random; it is an intrinsic stochastic process driven by a de- 
terministic external perturbation. This stochastic jumping leads to a much stronger noise 
in the power spectra. On the other hand, it is necessary that trajectories are within the 
vicinity of the separatrix of two different orbits so that the stochastic magnetization jump 
can take place under small perturbations. Since the separatrix has an infinite period, the 
trajectory close to it must have a nearly zero frequency as seen in Fig.BJd). 

The synchronization-like chaotic power spectra imply that quasi-periodicity can be ob- 
served in chaos if the trajectory is not near the separatrix. It is the stochastic jump of the 
trajectory in the vicinity of the separatrix that reduces the periodicity. The reduction of 
the periodicity due to stochastic jumps could also occur when the thermal fluctuations exist 
because the thermal fluctuations are stochastic, known as a Wiener process. In general, 

n 

there exists a bifurcation gap |24y where the thermal fluctuations limit the observation of 2 n 
orbits in the bifurcation diagram to a finite number n < n . The stronger the fluctuations 
are, the smaller n is. Consequently, the thermal fluctuations make the Lyapunov spectra 
smoother in Fig. 2 and chaos windows invisible in Fig. 3. We emphasize that although 
the thermal fluctuations lead to noise in the power spectra and the possible random jumps 
between two orbits, it is not associated with the bifurcation and thus there are no universal 
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Figure 4: Trajectories and their corresponding M z -component power spectra of two types of 
chaos. Upper panels: a,j = 262.5(Oe); the largest Lyapunov exponent is 0.86. Lower panels: 
dj = 261.0(Oe); the largest Lyapunov exponent is 1.56. All parameters are same as in Fig. 2. 



Feigenbaum numbers. Experimentally, one can distinguish chaotic dynamics generated by 
the deterministic perturbation (AC current) and by the thermal fluctuations. For exam- 
ple, when the parameter, e.g. DC current in our study, is varied, the current driven power 
spectrum should evolve from non-chaotic dynamics to chaos and back to non-chaos. 

In conclusion, by using both analytical and numerical approaches, we have shown the 
route to chaos and that the oscillator belongs to the class of dissipative quadratic non-linear 
systems. Our findings demonstrate that the current-driven magnetization dynamics is much 
richer than previously studied steady-state motion and synchronization. The magnetization 
oscillator whose dynamics can be measured by experiments provides a model system to 
verify chaos theories in a general non-linear system. 

This work was supported by DOE(DE-FG02-06ER46307). 
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